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CN ■ Abstract 
O 

^ ■ We have implemented a code for Z' + n jets production in ALPGEN, with 

Z' decays into several final states, including £ + £~ and tt. The MLM prescription 
' is used for matching the matrix element with the parton shower, including in 

this way the leading soft and collinear corrections. In order to demonstrate its 
; capabilities, we perform a combined analysis of Z' — > tt and Z' — > ttj production 

^ ' for a heavy leptophobic gauge boson. It is found that the effect of the extra jet 

(U | cannot only be accounted for by a if factor multiplying the leading-order cross 

section. In fact, the combined analysis for Z' — ► tt and Z' — > ttj presented 
improves the statistical significance of the signal by 25% (8.55c versus 6.77a for 
a Z' mass of 1 TeV), compared with the results of an inclusive analysis carried 
Q\ \ out on the same sample of tt + tij events. 
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1 Introduction 
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Searching for neutral vector resonances is one important task in the Large Hadron 

- T— I | 

^ ■ Collider (LHC) programme, being the Drell-Yan process pp — ► Z' — > £ + £~ , i = e, /i, the 

preferred one in these studies [1-4]. This signal with charged leptons in the final state 
has smaller backgrounds than those with only hadrons, but it requires that the new 
boson couples to both quarks and leptons, being this process suppressed when either 
of these types of couplings vanishes. The tt channel is an alternative if the couplings 
to leptons are the ones which are negligible. Among the hadronic final states, tt is 
very interesting by itself because, being the top very heavy, it allows for a relatively 
easy identification and reconstruction, and for this reason its backgrounds are relatively 
smaller. Even more, Z' — > tt is not only an alternative to the Drell-Yan process when 
the new resonance is leptophobic, but it is always complementary to determine the 
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model because it involves a different combination of couplings, and it also allows for 
asymmetry measurements in the semileptonic decay [5]. 

New gauge bosons are predicted in many of the best motivated Standard Model 
(SM) extensions. For instance, parity restoration, which can be at the TeV scale, 
requires extending the SM gauge symmetry to SU(2) L x SU(2) fi x XJ(1) B _ L , with new 
neutral and charged gauge bosons at this scale [6]. Grand unified models always predict 
new gauge bosons, and some of them may survive at lower energies [7]. In general, 
models addressing the hierarchy problem which are not based on supersymmetry, such 
as Little Higgs models [8] or models with large extra dimensions [9], also introduce (a 
plethora of) new vector bosons, with some of them at the LHC reach. Many of them 
couple to quarks and leptons but, as already stressed, even in this case (and obviously 
in the leptophobic limit) it is important to perform a detailed resonance search in the 
ti channel. 

With this purpose in mind, we have extended ALPGEN [10] with a generator for Z' 
production, including real emissions matched with leading logarithmic (LL) corrections 
(both real and virtual ones). This is presented in the next section, and applied to the 
search of vector resonances decaying into ti in the following ones. For the sake of 
illustration, we concentrate on the case of a leptophobic Z x , which also decays into 
new heavy neutrinos Z' x — > NN if m N < Mz> /2 [11]. The model is described in 
section [3], where we also give details of the simulation of top pair production mediated 
by neutral gauge bosons. Then, in section [4] we show the relevance of using a generator 
including these contributions: a combined analysis of Z' x — > ti and Z' x —> tij improves 
the LHC sensitivity to neutral vector resonances, raising the statistical significance 
of the signal by 25% compared with an inclusive analysis. Precise simulations of Z' 
production in different channels are not only essential for discovery, but to determinine 
the model by measuring all the Z' couplings with a precision as high as possible. The 
last section is devoted to our conclusions. 

2 The Z' generator 

In this section we briefly describe the main features of the new code. It evaluates the 
matrix elements for Z' + n jets production, with Z' decaying into several final states 
including and ti. Actually "Z' + n jets" stands for the sum of the three possible 
intermediate vector bosons, namely Z, 7* and Z'; their interferences as well as their 
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widths and all spin correlations are taken into account Detailed information can be 
found in the README file available at http : //mlm Thome . cern . ch/mlm/alpgen7] The 
SM parameters are controlled, as it happens for all the other processes implemented 
in ALPGEN [10], by the variable iewopt, whose default value is iewopt = 3. The 
possible final state is selected with the parameter ifs (e.g. if s = — > e + e~, if s = 
2 — > tt), and its left-handed (LH) and right-handed (RH) couplings to Z' with glzpe, 
grzpe and glzptop, grzptop respectively. In addition, the new gauge boson (LH 
and RH) couplings to the initial quarks, and its mass and width can be arbitrarily de- 
fined through the variables glzpup , grzpup, glzpdown, grzpdown, glzpc, grzpc, 
glzps, grzps, and masszp and zpwid, respectively, the lagrangian being 

£ = -i masszp (Z') 2 + ^ glzpup uj^l-^u Z'^+^ grzpup uj^l+j^u Z'^+. . . , (1) 

and zpwid the width of the Z' boson. Finally, the MLM matching prescription [12,13] 
can be applied by taking the parameter ickkw = 1. In Tabled] we gather the parton 
subprocesses evaluated for any given jet multiplicity. Since our aim is studying QCD 



Subprocesses 


uu -> // 


dd^ff 


gu -> ffu 


gd -> ffd 


ug -> ffu 


dg - ffd 


gg -> ffuu 


gg -> ffdd 



Table 1: Parton subprocesses (plus their charge conjugate) taken into account in the 
computation. Quarks u and d represent generic quarks of type up and down, respec- 
tively, and // stands for the fermion pair selected in the Z' decay through the variable 
ifs. If the jet multiplicity exceeds the number of light quarks in the final state, final 
state gluons are added up to reach the desired multiplicity. 



corrections to the production of a single s-channel Z' boson, rather than the production 
of two Z' bosons, we neglect contributions with two light quark pairs, because in our 
approach the former excludes the latter. 

In order to illustrate the capabilities of this code we will evaluate top pair production 
for the leptophobic model described in the next section. The full analysis is presented 
in section @J As it can readily be seen from Fig. [H taking M z > x = 1 TeV with a total 
width Y z> = 6.9 GeV, a proper calculation of the real radiation contributions to Z'j is 

1 The light jet(s) can also result from radiation off the Z' decay products, as for example in qq — > 
Z' — > tig, but for simplicity we still denote these processes as Z' + n jets. Obviously, by Z' — ► ttj we 
only refer to events in which the extra jet is radiated off the Z' decay products. The radiation from 
top decay products is not included at the matrix element level. 
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Figure 1: tt versus ttj invariant mass distributions at the generator level, for Z' x j events 
originated in qq (left) and gq (right) collisions at LHC. The number of points in the 
plots (8172 for qq and 1332 for gq) corresponds to the expected number of events at 
the LHC for M z > = 1 TeV with a luminosity of 30 fb _1 . 



compulsory to account for the relatively large number of Z' — > ttj events produced in 
qq collisions J (In both cases the SM contributions from Z and 7* are turned off, but 
they are included in the analyses below.) In particular, as it can clearly be seen in the 
left plot, 37% of the qq — > ttj events correspond to a ttj resonance, being the extra 
jet from final state radiation (FSR) from one top quark. The rest corresponds to a tt 
resonance with the jet from initial state radiation (ISR). 



3 Description of the model and the simulation 

As explained in the previous section, the code allows for arbitrary Z' couplings to 
fermions (models). For illustration we choose a model with vanishing couplings to 
ordinary leptons, but non-zero couplings to SM quarks (in particular to the top) and to 
new heavy neutrinos. Being the new boson leptophobic, the lower bound on its mass is 
rather weak, and sizeable signals are already possible at Tevatron [11]. There are many 
models with extra leptophobic gauge bosons, originally studied to interpret the initial 
disagreement between the LEP data onZ^ bb, cc [15] and the SM predictions [16,17]. 
We will restrict ourselves to an E 6 based model [18]. The neutral gauge interactions 
are described by the Lagrangian [19] 



T 3 gW£+y/lYg Y B» + Q , g'z'» 



V>, (2) 



2 Notice that in the SM, when usual isolation and transverse momentum cuts are applied, the largest 
contribution to Zj comes from gq collisions [14]. 
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where a sum over the Weyl fermions of the fundamental E$ representation 27 and 
the three families must be understood. Y is the SM hypercharge properly normalised, 
and the extra Z' x charges Q' correspond to the only leptophobic combination within 
E 6 [17,20], Q' = 3/y/W (Y v + Y/3) , with Y v the extra U(l) defined by flux breaking [21]. 
For the SM (LH) fields 

2Q U = 2Qd = Qu c = ~^Qd c = ~~ 5 

Q'v = Q'e= Q'e° = , (3) 

reading the code charges for ilep = 2 

glzptop = glzpup = glzpc 
grzptop = grzpup = grzpc 



2^/6 

VE' 



glzpdown = glzps = grzpdown = grzps = -= . (4) 

2v6 

A detailed discussion of the phenomenology of this SM extension can be found in 
Ref. [17], where a nearly-leptophobic model with Q' ~ Y v + 0.29 Y is studied together 
with several other alternatives. In our phenomenological study we will assume for 
simplicity that the extra vector-like lepton doublets and quark singlets of charge —1/3 
are heavier than M Z ' x /2, as they are the heavy neutrinosjf] Possible supersymmetric 
partners are taken to be heavier as well. Otherwise, the total Z' x width would be 
larger, decreasing the cross sections into SM final states. The total leading order (LO) 
Z' x production cross section at LHC is plotted in Fig. [2] as a function of Mz' , together 
with the maximum (i.e. when Z' x decays only to SM fermions) cross section for tt 
final states. The coupling constant of the new U(l)' has been fixed for reference as 
g' = a/5/3 gy = a/5/3 g sw/cw, and cross sections are calculated using CTEQ5L 
parton distribution functions [23]. 

For the simulations in the next section we also take a Z' x mass of 1 TeV, above 
the Tevatron exclusion limits for this model [24]. The signal is generated with Monte 
Carlo statistics of 300 fb _1 and rescaled to 30 fb _1 . SM backgrounds include tinj (with 
nj standing for n jets at the parton level), single top, W/Znj, W/Ztinj, W/Zbbnj, 
W/Zccnj, diboson and triboson production. They are generated with a Monte Carlo 
statistics of 30 fbr 1 using ALPGEN, taking m t = 175 GeV, M H = 115 GeV. (The 
complete list of processes and numbers of events generated can be found in Ref. [25].) 



3 If Z' x can decay into heavy neutrino pairs NN (with Z' x charge Q' N = — 3/2V6), this could be 
also observed in multi- lepton channels [11,22]. 
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Events are interfaced to Pythia 6.4 [26] to add soft ISR and FSR and pile-up, and 
perform hadronisation. The MLM prescription is also used to perform the matching 
for the backgrounds, with default values for the matching parameters. A detailed 
investigation of the uncertainties related to the matching procedure with ALPGEN 
for the main background processes has been presented in Refs. [13,27]. In order to 
simulate a real detector environment we use the fast simulation AcerDET [28] which 
is a generic LHC detector simulation, neither of ATLAS nor of CMS, with standard 
settings. In particular, the lepton isolation criteria require a separation AR > 0.4 
from other clusters and a maximum energy deposition HE? = 10 GeV in a cone of 
AR = 0.2 around the reconstructed electron or muon. Jets are reconstructed using 
a cone algorithm with AR = 0.4. In this analysis we only focus on central jets with 
pseudo-rapidity \r]\ < 2.5. For central jets, a simple b tagging is performed with 
probabilities of 60% for b jets, 10% for charm and 1% for light jets. 

Our estimate of the signal relies on the LO approximation. In the absence of a 
proper next-to-leading order (NLO) calculation, we estimate the impact of hard NLO 
corrections using MCFM v5.6 NLO code [29]. We have rescaled Z and W masses 
(and accordingly the Gp coupling constant) in such a way that the Z mass is pushed 
to 500 GeV and 1 TeV, while keeping the SM couplings, in order to simulate the 
sequential Z' SM . We selected the bb final state among the possible ones, choosing a 
window Mz> — 7T Z > < m h j. < M z > + 7T Z ' ■ As a result the i^-factor, defined 
as a(NLO)/cr(LO) with renormalization and factorization scales fixed to Mz' SM , turns 
out to be of the order of 20% for both M z > = 500 GeV and M z > = 1 TeV. 



6 



4 Z' —> tt and Z' —> ttj observation at LHC 



The presence of a Z' resonance is detected as a bump in the tt invariant mass spec- 
trum, which can eventually be normalised from real data with a good precision. For 
this analysis we restrict ourselves to the tt semileptonic decay channel, in which the 
kinematics can be fully reconstructed. The pre-selection criteria are: (i) one charged 
lepton with pr > 30 GeV; (ii) two 6-tagged jets with pt > 20 GeV; (iii) at least two 
light jets with px > 20 GeV; (iv) total transverse energy H T > 750 GeV. The latter 
cut is implemented in order to reduce the QCD ttnj background, as well as to remove 
from the signal processes the SM component mediated by Z and 7*. Note that these 
cuts can be optimised to enhance the statistical significance of the signal but, on the 
other hand, reducing the background too much makes its normalisation from real data 
difficult. The number of events fulfilling these requirements are collected in Table [2] 
(left) for the signal and main backgrounds (the rest of backgrounds are not explicitly 
shown but we keep them in the calculations). On the right panel we show the decom- 
position of the ttnj background in n = 0, . . . , 5 multiplicities. The dominance of the 
n— 1,2,3 subsamples results from the H T cut, which has a much higher rejection for 
lower jet multiplicities. 













Pre. 


acc. Ht 




Pre. 




Pre. 


tiO] 


7512 


2.7% 




493.2 




1213.9 


til] 


16665 


9.1% 


ttnj 


64688 


ttbb 


657 


tt2j 


16392 


22.7% 


tw 


2425 


Wnj 


915 


mj 


10299 


43.6% 


tj 


445 


Wbbnj 


2202 


tUj 


4580 


66.6% 










tibj 


1601 


84.6% 



Table 2: Left: number of signal and background events at the pre-selection level. Right: 
decomposition of the ttnj background in subsamples with n jets at the partonic level 
and acceptance of the H T cut. The luminosity is 30 fb _1 . 

The first step to search for the Z' signal is the reconstruction of the top quark 
pair. These are found by choosing the best pairing between b jets and reconstructed 
W bosons: 

1. The hadronic W is obtained with the two jets (among the three ones with largest 
Pt) having an invariant mass closest to Mw- 

2. The leptonic W is obtained from the charged lepton and the missing energy, 
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Figure 3: Reconstructed top quark masses. The luminosity is 30 fb 



identifying (p v )r = tfr , requiring (pe+p u ) 2 = and solving for the longitudinal 
component of the neutrino momentum. If no real solution exists, the neutrino 
transverse momentum is decreased in steps of 1% and the procedure is repeated. 
If no solution is still found after 100 iterations, the discriminant of the quadratic 
equation is set to zero. Both solutions for the neutrino momentum are kept, and 
the one giving best reconstructed masses is selected. 

3. The two top quarks are each reconstructed with one of the W bosons and one 
of the b jets, and are labelled as 'hadronic' and 'leptonic' corresponding to the 
hadronic and leptonic W, respectively. 

4. The combination minimising 

(m^ d -M w ) 2 {m^-M w f (m^ - m t f (mg - m t f 

2 ' 2 ' 2 2 V" 3 / 

is selected, with a w = 10 GeV, a t = 14 GeV [2]. 

We present the reconstructed mass distributions at the pre-selection level in Fig. [31 
With the top quark pair identified, one can consider several variables to discriminate 
the signal from the background. The most interesting ones are the ti, tij invariant 
masses and the transverse momentum of the ti pair. They are presented in Fig. [4] for 
the two signal subsamples with zero (Oj) and one (lj) jet at the partonic level and 
for the SM background. For completeness, we also show the transverse momentum 
distribution for the top quark decaying leptonically (for the hadronic top quark it is 
similar), which is not used to enhance the signal significance The most remarkable 
features, which guide our subsequent analysis, are: 



4 Although at pre-selection this variable seems to have a good discriminating power for Z' (Oj) 
events, after event selection criteria based on and other variables the signal and background 
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Figure 4: Kinematic distributions at pre-selection. Up, left: tt invariant mass; up, 
right: ttj invariant mass (requiring at least three light jets); down, left: transverse 
momentum of the tt pair; down, right: transverse momentum of the top quark decaying 
leptonically (not used for event selection). The luminosity is 30 fb _1 . 



1. The m t i distribution has a large off-peak component from the Z' (lj) sample. If 
the SM background cannot be predicted with very good accuracy (as it seems the 
case for tt plus several hard jets), and must be normalised from off-peak data, 
this tail will behave as a combinatorial background reducing the peak significance 
unless a specific analysis is carried out with a different reconstruction for these 
events. 

2. The Z' (lj) sample exhibits a good peak in the tij invariant mass distribution, 
although slightly shifter towards m t fj values larger than 1 TeV. 

3. We also observe that the tt transverse momentum is typically much smaller for the 

distributions become very similar. Nevertheless, this and other distributions, as for example the tt 
pair rapidity, can be implemented in a likelihood function to achieve a better discriminating power 
than with the cut-based analysis implemented here. Such optimisation is beyond the scope of the 
present work. 
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Z 1 (Oj) signal than for the background. Although tt pairs are typically produced 
with low transverse momentum, the requirement > 750 GeV has a higher 
suppression for lower hard jet multiplicities, as seen in Table El As a result of 
this cut, the distribution is shifted towards large values, as it can observed in 
Fig. El This higher jet multiplicity is also the reason for the worse reconstruction 
of the hadronic top in Fig. [31 
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Figure 5: tt transverse momentum distribution for the ttnj background at the pre- 
selection level with and without the Ht cut. The luminosity is 30 fb -1 . 



The best strategy to maximise the signal significance is to develop two analyses. The 
first one aims to reconstruct the it peak in the region of low transverse momenta 
of the tt pair setting pj, < 50 GeV, which eliminates a large fraction of background 
and the off-peak signal contribution which would otherwise constitute a combinatorial 
background. The second analysis will search for the ttj peak in the complementary 
region > 50 GeV. We present these two analyses in turn. Then, for comparison, we 
will perform an inclusive analysis without separate reconstructions. 

4.1 Analysis I 

As selection criteria we impose some loose quality cuts on reconstructed top quark 
masses and, more importantly, small transverse momentum for the tt pair, 

125 GeV < m^ ad ,m' ep < 225 GeV, 

V% < 50 GeV . (6) 

The number of signal and background events with these cuts can be read in Table d 
We observe that the p% cut significantly reduces the ttnj background and removes 90% 
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Sel. Peak Sel. Peak 

Z' x (Oj) 334.5 299.2 Z' x (lj) 149.0 114.4 



ttnj 


6581 


1652 


tibb 


15 


3 


tw 


85 


17 


Wnj 


39 


13 


tj 


9 


4 


Wbbnj 


15 


4 



Table 3: Number of signal and background events at the selection level (analysis I) and 
at the tt resonance peak. The luminosity is 30 fb _1 . 



500 




Figure 6: tt invariant mass for the SM background and the background plus the Z' 
signal at the selection level (analysis I). The luminosity is 30 fb -1 . 



of the Z' signal with n — 1, in particular the off-peak contribution. The presence of 
the Z' resonance can be spotted with the analysis of the tt invariant mass distribution, 
presented in Fig. [6] for the background alone and the background plus the signal at the 
selection level. The number of events at the peak 

900 GeV < m t t < 1100 GeV (7) 

can be read in Table [31 We will assume that we normalise the background by off peak 
measurements, obtaining a scaling factor k = 1.05 determined from the comparison of 
the distributions for signal plus background and background alone. Then, the statistical 
significance of the peak is S' / y/B', where 

B' = kB S' = (S + B) - B' (8) 

and S, B are the true numbers of signal and background events. For the peak region 
in Eq. ([7]), the excess of events amounts to 7.78<r. 
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4.2 Analysis II 



The main motivation for this analysis is the fact that many of the Z' (lj) events do not 
exhibit a peak in the it invariant mass distribution and would fall off the peak region 
in Eq. ((Jj). This is seen clearly in Fig. 0] (up, left). When the invariant mass of the it 
pair plus the hardest additional jet is considered (requiring in this case a minimum of 
three jets) the distribution exhibits a clear peak, although slightly displaced, as shown 
on the upper right panel of that figure. We then concentrate on the complementary 
event sample for this analysis, requiring at least three light jets and setting the cuts 

125 GeV < m^ ad ,m' ep < 225 GeV, 
■p% > 50 GeV , 

mm(A J R Mmd ,A J R Mep )<1.6 (9) 

as selection criteria. The last one is implemented to reduce the background, since the 
signal events peaking at m t ij ~ Mw are produced by FSR in Z' — > it. The number 
of events after these cuts are given in Table 0J The itj distribution for the signal plus 



Sel. Peak Sel. Peak 

Z' x (Oj) 52.2 41.7 Z' x (lj) 412.8 293.5 
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Table 4: Number of signal and background events at the selection level (analysis II) 
and at the ttj resonance peak. The luminosity is 30 fb _1 . 

background and background alone is presented in Fig. The number of events at the 
peak 

900 GeV < m t ij < 1200 GeV (10) 

can be read in Table [H The background scaling factor in this case is k = 1.022, and the 
excess of events over the (normalised) SM background expectation amounts to 3.56<r. 

4.3 Inclusive analysis 

In order to see the advantage of separate, dedicated analyses for Z' — > ti, Z' — > ttj, we 
also perform the inclusive analysis for the two Z' signal components, searching for a 
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Figure 7: ttj invariant mass for the SM background and the background plus the Z' 
signal at the selection level (analysis II). The luminosity is 30 fb -1 . 



peak in the m t t distribution. For event selection we only require a good reconstruction 
of the top quark pair, 

125 GeV < m t had , m' ep < 225 GeV , (11) 

and drop the other kinematic cuts. The number of signal and background events can 
be read in Table [5l 

Sel. Peak Sel. Peak 

Z' x (Oj) 400.1 332.6 Z' x (lj) 727.1 361.0 

ttnj 35203 4196 tibb 163 9 

tW 567 40 Wnj 170 20 

tj 109 12 Wbbnj 263 13 

Table 5: Number of signal and background events at the selection level (inclusive 
analysis) and at the tt resonance peak. The luminosity is 30 fb -1 . 

The tt distribution for the signal plus background and background alone is presented 
in Fig. El The number of events at the peak 

900 GeV < mu < 1100 GeV (12) 

is also given in Table [51 As expected, in this case the background scaling factor is 
larger than for the other analyses, k = 1.055, and the peak significance is of 6.77<r. 
This number is 25% smaller than the combined significances of the other two analyses, 
8.55a. 



13 



2000 




2000 



">,1 



Figure 8: tt invariant mass for the SM background and the background plus the Z' 
signal at the selection level (inclusive analysis). The luminosity is 30 fb _1 . 



5 Conclusions 

We provide a new code implemented in ALPGEN for neutral vector resonance pro- 
duction, including the higher order corrections from real emission and the virtual LL 
contributions. In order to illustrate its features we have generated ti and tij events by 
the exchange of a leptophobic gauge boson based on E 6 at LHC. The analysis of the 
Z' x — > tt and Z' x — > tij samples show that such a program is necessary to account for 
~ 23% of the tiX events. Indeed, for these events the resonance is not found in the 
tt invariant mass distribution but on the tij invariant mass. Therefore, their presence 
cannot be accounted for by a K factor multiplying the LO Z' x —> ti cross section but 
has to be properly simulated at the generator level. 

We have shown that the LHC discovery potential for neutral vector resonances 
in the tiX final state benefits from a separate analysis for ti and tij. For the case 
examined (with a Z' x mass of 1 TeV), the enhancement over an inclusive search for ti 
resonances is of a 25% in the statistical significance. This improvement is expected to 
be maintained for larger masses, and of course for other Z' models than the leptophobic 
Z' x one used in the simulations. For other ti resonances the trend is expected to be the 
same as well. 

Finally, the inclusion of corrections to Z' — > ti not only translates into a better 
discovery potential but also into a proper description of the kinematical distributions. 
A code implementing them is necessary to establish the model once a new vector 
resonance is discovered, discriminating among models by a more precise measurement 
of the angular distributions [30] and the different Z' couplings [8,19]. 
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